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THRUSTING TRAJECTORY MINIMIZATION PROGRAM 
FOR ORBITAL TRANSFER MANEUVERS 


By Lawrence H. Hoffman, Richard N. Green, 
and George R. Young 
Langley Research Center 

SUMMARY 

A computer program has been designed which determines the minimum -burn -time, 
thrusting, transfer trajectory between two Keplerian orbits. The minimization equations 
are formulated with constant Lagrange multipliers and solved numerically with the 
Newton -Raphs on method. The solution obtained in this paper is not truly optimum 
because the control vector has been restricted to constant values (i.e., to be optimum the 
control vector should be a function of time). 

The equations of motion for the transfer trajectory are those of a spacecraft maneu- 
vering with constant thrust and mass -flow rate in the neighborhood of a single body. The 
thrust vector is allowed to rotate in a plane with a constant pitch rate. The transfer 
trajectory is characterized by six control parameters and the final orbit is determined or 
partially determined by the desired target parameters. The program is capable of 
varying from one to six control parameters to find the desired target parameters which 
are chosen from a large set. If the number of target parameters is less than the number 
of control parameters, the fuel required for the maneuver is minimized. To conserve 
computer time the equations of motion are integrated by a truncated power series in time. 
The use of the program is illustrated with three sample computer cases. 

INTRODUCTION 

In order to study completely an interplanetary, orbital mission it is necessary to 
define the maneuver targeting logics that are required at each of the guidance junctures. 

In the case of the Viking mission it is felt that adequate interplanetary targeting studies 
can be accomplished with the Mark IV Error Propagation Program (ref. 1) or the 
Simulated Trajectories Error Analysis Program (ref. 2). Both of these programs include 
targeting options to determine the midcourse velocity corrections such that the trajectory 
constraints are satisfied at the planet. In addition, each program can be used to perform 
an error analysis of the interplanetary phase of the mission. However, since neither 
program has a finite burn maneuver capability, they will not solve the targeting problem 
for the Mars orbit insertion maneuver, the orbit trim maneuvers, or the deorbit maneu- 
vers. Various programs exist for the near planet targeting analysis required for this 



mission (e.g., refs. 3, 4, and 5); however, for various reasons the existing programs were 
unsuitable for this analysis. References 3, 4, and 5 are all for a coplaner transfer which 
is not general enough for the Viking mission. Many other existing programs are impul- 
sive and therefore one cannot estimate fuel requirements accurately with them. For 
these reasons a computer program VITAP (Viking Targeting Analysis Program) was 
developed which solves the near Mars phase of the targeting analysis. 

One of the areas of investigation in Project Viking is that of determining the 
minimum -fuel, thrusting, transfer trajectory between two Keplerian orbits. This pro- 
blem exists for the Mars orbit insertion maneuver, the orbital trims, and for the deorbit 
maneuver. The equations of motion for the transfer trajectory are those of a spacecraft 
maneuvering with constant magnitude thrust and mass -flow rate in the neighborhood of a 
single body. To increase the flexibility of the program the thrust vector is allowed to 
pitch at a constant rate; however, for Project Viking, the pitch rate is zero. Therefore, 
the problem considered here is as follows: Given an initial Keplerian orbit and the 
characteristics of the engine, determine a set of control parameters which define the 
thrusting maneuver such that the required fuel is minimized and the resulting orbit 
satisfies a number of constraints. 

A description of the mathematics used in VITAP is presented. It consists of solving 
a finite -dimensional minimization problem with equality constraints. The minimization 
equations are formulated with constant Lagrange multipliers and solved numerically with 
the Newton-Raphson method. Since these equations are very complicated, the first and 
second partial derivatives necessary for the solution are computed numerically. A dis- 
cussion of this procedure is also included. A detailed description of the assumed model 
of the thrusting spacecraft is given together with the associated equations of motion. To 
conserve machine time these equations are integrated by a truncated power series in 
time. The power series solution to the equations of motion is completely developed in 
the appendix. A description of the computer program is given with a discussion of the 
program input and output. The use of program VITAP is illustrated with three sample 
computer cases. 


SYMBOLS 

A submatrix defined in equation (8) 

a semimajor axis, kilometers 

a* specific value of a, kilometers 

a(t) magnitude of thrust acceleration at time t, kilometers/second2 
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B submatrix defined in equation (8) 

B vector from center of planet to incoming hyperbolic asymptote, kilometers 

B-R component of B in R direction, kilometers (see sketch 3) 

B*T component of B in T direction, kilometers (see sketch 3) 

e eccentricity 

F augmented function (see eq. (1)) 

f(5) function of a 

G lighting angle at landing point, degrees (see sketch 2) 

g partial derivatives of F with respect to a (see eq. (4)) 

i inclination, degrees 

i* specific value of i, degrees 

K step size control variable (see eqs. (11)) 

m mass, kilograms 

N number of trajectory integration increments 

R error variable (see eq. (12)) 

R unit vector perpendicular to planet equator (see sketch 3) 

r radius from center of planet, kilometers 

r^ radius of apoapsis, kilometers 

rp radius of periapsis, kilometers 

S unit vector parallel to incoming hyperbolic asymptote (see sketch 3) 
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T 

t 

tb 

Vco 

% 

x,y,z 

a,/3,6 

a 

“1 

AV 


unit vector in planet equator perpendicular to S (see sketch 3) 
time, seconds 

time duration of thrusting maneuver, seconds 
hyperbolic excess velocity, kilometers/second 
six-dimensional state of spacecraft at end of maneuver 
rectangular Cartesian coordinates, kilometers 
angles defining direction of thrust, degrees (see sketch 1) 
vector of control variables 
first guess for a 

integral of acceleration due to thrust, kilometers/second (see eq. (18)) 


Aa = a - aj 


AaijAd!^ components of vector A«* 
a 5* maximum allowable step size of a 


AX - Xj 

6a infinitesimal variation of a 

e convergence criteria (see eq. (12)) 

6 angle between landing point and periapsis, degrees (see sketch 2) 

6 thrusting pitch rate, degrees/second (see sketch 1) 

X vector of constant Lagrange multipliers 

Xj first guess for X 
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jLi gravitational constant of planet, kilometers V second^ 

V true anomaly, degrees 

4> latitude of landing point, degrees 

rl/ trajectory constraints 

O longitude of ascending node, degrees 

Cl) argument of periapsis, degrees 

Subscripts: 

f final value 

i,j,Z indices 

m number of constraints ^ 

mXm matrix with m rows and m columns 

min minimum 

n number of control parameters 

o initial conditions 

1,2,. . . first, second, . . . 

Superscripts: 

T matrix transpose 

-1 matrix inverse 

’ modified parameter 

Dot over a symbol indicates differentiation with respect to time. 
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ANALYSIS 


Finite -Dimensional Minimization Problem with Equality Constraints 

Consider the finite -dimensional minimization problem in which equality constraints 
have been imposed. It is desired to find a set of controls . .,Q!n) 

function • *>“ 11 ) is minimized subject to the m constraint conditions 

^lia) = 0 
i/ZgCc?) = 0 


’^m^^ = 0 

where m g n. If f and are sufficiently smooth and linearly independent, the 
solution must satisfy the conditions 


da 
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8X 
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8F 
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9«n. 


8X]^ 
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where F is defined as 
Z=m 

F =f + 

1=1 
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( 2 ) 
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and 




8F 

8X 






= 0 


(3) 


mj 

In addition, a sufficient condition for a local minimum (refs. 6 and 7) is given by 
i=n i=n 

6q'?6q'j > 0 

Sofi dai J 


i=i i=l 


'1 


where 


da = "a - a 


mm 


and is constrained to satisfy 




9a!j 

aofn 




8ofn 


6oti 


6 Of 


n 


= 0 


Newton-Raphson Technique 

The necessary conditions for a minimum (eqs. (2) and (3)) can be solved by the 
Newton-Raphson technique. Let 


(4) 


Then the necessary conditions can be written as 
gi(5,l) = 0 
^lCS) = 0 


(i = l,2,...,n) (5) 

(i = l,2,...,m) (6) 
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Expanding equations (5) and (6) in a truncated Taylor series about the point 

(where denotes the first guess at the control vector and Xj denotes the first guess 

at the Lagrange multipliers) yields 


gi(3,x) = gi(aiAi) + - «l) + — -\) -0 (i - 1,2,. . .,n) 


hi'^) = H^l) 



{I - 1,2,. . .,m) 


Since 

_ 3 /8F(a,x)\ _ 

dXi dXi\ Q“i / 80 !i 

then the necessary conditions can be expressed as 
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where Aa s 3 - and Alt = X - are the corrections to improve the initial guess. 
The familiar Newton-Raphson iteration equations can now be written as 
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Inspection of equation (7) reveals that it can be written as 
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( 8 ) 


For the case where m = n, it is instructive to consider the inverse of the matrix in 
equation (8) using Schur's method for partitioned matrices (ref. 8), that is, 


Aotj 


A"^ - A-^b'^(bA-1b'^)'^BA"^ 

I 

; a"^b’^(ba'^b’^) ^ 
1 
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If m = n and all the proper inverses exist, the expression A“^ - (bA'^b"^) BA 

reduces to the null matrix [(^ and the corrections to the control parameters depend only 
on the errors in the constraints and not on the Lagrange multipliers associated with 
gj(o!,x). This is to be expected when n controls are varied to satisfy n constraints 
and no minimization takes place. In addition, when m = n the expression 

,-l 


-1 
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reduces to B~ such that 
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( 10 ) 
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which are the Newton-Raphson iteration equations for solving n equations in 
n unknowns. 

Experience has shown that the range of convergence for the Newton-Raphson tech- 
nique can be considerably extended by placing constraints on the allowable size of the step 
to be taken. If the step size of the controls is constrained to be less than Aa , then a 
modified step size is given by 


where 


Affj = K"^ AcKj 
AX'^ = K"^ AX; 


K = max 


|A«i| |A«2| 

^ Affi Ao-n 


• J 


^“n / 


(i = 1,2,.. .,n) 
(1 = 1,2,. . .,m) 


( 11 ) 


Such a procedure will greatly increase the range of convergence for equation (9). 

The Newton-Raphson technique is considered complete when the errors in the con- 
straints and the first partial derivatives of the augmented function are small, that is, when 
the control vector converges and ceases to change. Symbolically, the solution is complete 
when 


R < e 


( 12 ) 


where 

R = |gj| + . . . + |gn| + + • • • + 

and 6 is a small number. 


Numerical Differencing Formulas 

The Newton-Raphson iteration technique used to solve the finite -dimensional mini- 
mization problem requires expressions for a number of partial derivatives. From equa- 

8g. 

tion (7) it can be seen that the partial derivatives and —1 are needed where 

6«j 

80!i 90 !^ 9 Q!. Zj ^ 80!. da. 

J J 1=1 1 J 

Assuming that f(afj,Qf 2 ,. . -j^n) 1*® differentiated directly, it remains to obtain the 

first and second derivatives of the constraint functions with respect to the con- 
trols In general • v<^n) 1® ^ ® and does not lend 
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itself to direct differentiation. However, it was found that numerical partial derivatives 
are quite adequate for the Newton-Raphson iteration equations. Letting Aai denote a 
small increment in o!i, the first partial derivatives can be represented by (ref. 9, p. 136) 


IS = • -’“n) - h («1- • • -«n)~ 

J J 


(13) 


where I = 1,2,. . .,m and j = 1,2,. . .,n. Similarly, the second partial derivatives can be 
represented by 




• v«i+A0!i,. . .,Q!j+Aa!j,. . .,Q!n) 


Bai 9cKj 4Ao'i Aojj 

+ {oil,- • • -»“n) (14) 

where I = 1,2,. . .,m and i,j = 1,2,. . .,n. When j is replaced by i, equation (14) 
reduces to 


1 




H (“!’* • • ->“n) - (<^l»- • v“n) + h (“l>- • vai-Aofi,. . .,«n) 


which reduces the number of times that is evaluated. It can be shown that 

must be evaluated 2mn times to obtain the first partial derivatives and f ^ f ^^'5 77 + 1 

\2I (n - 2). , 

times to obtain the second partial derivatives. 


APPLICATIONS 


Equations of Motion 

The engine used to perform these maneuvers is assumed to have a constant thrust 
and a constant mass-flow rate m. If the mass of the spacecraft at the start of the 
maneuver is mQ, then the mass at time t is given by 

m(t) = mQ + fht (15) 
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and the magnitude of the acceleration by 


a(t) = 


Thrust 
mo nat 


The direction of the acceleration vector is defined by sketch 1, where the X,Y,Z triad is 
the inertial Cartesian coordinate system alined such that the Z-axis is perpendicular to 
the Mars equator and the X-axis is at the Mars vernal equinox. This system of coordi- 
nates is normally referred to as the areocentric system. From sketch 1, it can be seen 


Z 



that the locus of the thrust vector is in the plane defined by the angles « and /3. At the 
start of the maneuver (t = 0) the thrust -vector direction is defined by the angle 6 and 
thereafter is allowed to rotate in the plane at a constant rate 0 until the maneuver 

is terminated at t = tt,. E>uring the maneuver the thrust vector is defined by the angle 
6+0t and the direction cosines are 


cos (6+0t) cos Q? - sin (6+0t) sin a cos )3 
cos (6+Ot) sin a + sin (6+0t) cos a cos j3 
sin (6+0t) sin j3 

The assumed trajectory model is two body motion plus an acceleration due to thrust and 
is defined by the equations of motion 
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X = + a(t) [cos (6+0t) cos a - sin (640t) sin a cos 

( 16 a) 

y = — ^ + a(t) [cos (6+0t) sin a + sin (6+0t) cos a cos 

( 16 b) 

z = + a(t) [sin ( 5 + 0 t) sin 

( 16 c) 


where 

a(t) = J3>£HfL 

mQ + mt 

r =(x^ + + z^) ^ 

and jLi is the gravitational constant. The initial conditions for the equations of motion 
are derived from the knowledge of the initial orbit plus a control variable which 
denotes the true anomaly on the orbit at the start of the maneuver. In other words, the 
initial orbit is known, but the "best" position along this trajectory to perform the maneu- 
ver is unknown. Therefore, the best position Uq must be determined during the optimi- 
zation process. In fact there are six parameters that characterize the maneuver and 
must be determined, namely, a,^, 5 , 9 The initial direction of thrust is defined 
by O', 0 , 6 ; the pitch rate, by 6 ; the time duration of the burn by t^; and the position 
of ignition, by Vq. These six parameters constitute the set of control parameters 
^Qtj,Q!2,. . -,Q!n) finite -dimensional minimization problem. The function of these 

parameters to be minimized is simply 

• -.“n) = = t^j^ 

Since the mass -flow rate is constant, minimizing the square of the bum time is equivalent 

0 

to minimizing the required fuel. The function to be minimized is not restricted to tj^ 
and there are other continuously differentiable functions of the controls that could have 
been used. Once the six controls are determined they along with the initial orbit com- 
pletely define the final orbit. 


Trajectory Constraints 

The controls are to be determined such that the burn time is minimized subject to 
the requirement that the final orbit has certain characteristics. Therefore, the control 
parameters are required to satisfy certain constraints, that is 
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( 17 ) 




where m = 6. Since the maneuver is defined by six control parameters, in general no 
more than six constraint functions can be satisfied. For the burn time to be minimized, 
the number of constraints must be less than the number of control parameters. As an 
example it might be required that the final orbit have a semimajor axis a equal to a 
constant a* and an inclination i equal to i* . The constraint functions would then be 

= a(3^) - a* = 0 

^2=i(X() -i* =0 

where a and i are functions of the state Xj at the end of the burn and is a 
function of the controls through the equations of motion (eqs. (16)). Since the values of 
a and i are fixed in the final orbit, they are referred to as target parameters. The 
minimization process targets to the specified values of these parameters while minimizing 
the burn time. 

There are any number of target parameters from which to choose. An obvious set 
is the six Keplerian orbital elements a, e, i, co, and v. Others would include the 
radius of periapsis rp and the radius of apoapsis rg^. The period of the final orbit is 
a likely candidate. It might be required to establish the final orbit in such a manner that 
the spacecraft would be at a given position at some future time. Therefore, the latitude, 
longitude, and true anomaly of the spacecraft at a reference time should be included in 
this list of target parameters. 

Consider the requirement that the landing point on the surface be located at a pre- 
scribed lighting condition or at a given latitude. The geometry of this situation is shown 
in sketch 2. 
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The landing point is arbitrarily defined to be at an angular distance 9 from periapsis 
and within the plane of the final orbit. Since 9 is assumed constant, the landing point is 
a function of only the orientation angles i,o),S2 of the final orbit. Therefore, the 
thrusting maneuver could be constrained to an orientation that places the landing point at 
a prescribed latitude 4>- Similarly, the lighting angle G at the landing point could 
serve as a target parameter. 

Other useful target parameters are the h 3 ^erbolic impact plane param- 
eters B*T, B*R, Vco, S. These parameters are defined in sketch 3 and reference 10. 

The three unit vectors S,T,R form a coordinate system where S is parallel to the 
incoming asymptote of the spacecraft orbit, T lies in the Mars equatorial plane perpen- 
dicular to S, and R completes the triad. The miss vector B is in the R-T plane 
and represents the distance from the center of Mars to the incoming asymptote. It is 
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Sketch 3 


usually characterized by its two components B-T and B-R. This set of parameters 
B-T, B-R, S and the hyperbolic excess velocity completely specify the final 
hyperbolic orbit. They are frequently more useful than the standard orbital elements 
and, therefore, should be included in the list of target parameters. 

The choice of specific target parameters determines the constraint equations which 
enter the optimization process. As stated previously, no more than six constraints and at 
least one constraint must be imposed. The choice of the individual target parameters is 
important in that the ones mentioned are not all independent. As an example, suppose 
that the three target parameters a,e,rp were chosen. Since a and e dictate the 
radius of periapsis of the final orbit, it would be impossible to satisfy these three con- 
straints unless the third constraint happened to be consistent with that dictated by the first 
two. In general these three target parameters are dependent and would not constitute an 
acceptable set of constraints. Therefore, the specific target parameters should be chosen 
with care. 
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The areas in Project Viking which require optimum orbital transfer can now be 
analyzed. Given the orbital elements of the original orbit and an initial guess on the 
controls and the multipliers , the necessary equations for a minimum 

maneuver (eqs. (2) and (3)) subject to the desired constraints can be solved by the 
Newton -Raphson iteration technique (eq. (7)). 

COMPUTER PROGRAM VITAP 
Program Description 

The equations for an optimum, thrusting transfer between two Keplerian orbits and 
their associated Newton -Raphson solution have been incorporated into a computer pro- 
gram, VITAP. It consists of a main program and 11 subroutines which are written 
entirely in FORTRAN computer language for the Control Data 6600 computer system. 

The program resulted in a field length of SOOOOg. 

Various options have been included in VITAP which allow for a considerable amount 
of flexibility. The option is available to vary the six control parameters 
or to fix one or more of the controls at a constant value. For example, if the pitch rate 
6 is fixed at a constant value, the transfer maneuver is optimized with respect to the 
remaining five control parameters. This option allows various guidance laws to be 
considered. In addition a nvimber of different target parameters are available as con- 
straints on the optimum maneuver. In fact, as many as six constraints out of a set of 
twenty can be selected. Thus, quite a large number of combinations of control variables 
and target parameters are available. Program VITAP also operates in three modes. 

The first mode is the normal optimum transfer outlined previously. The second mode 
allows for the optimum transfer with the additional constraint that the inclination of the 
final orbit be between an upper and lower bound. The third mode targets backwards. 

This mode considers the problem of finding the best orbit from which to establish a given 
orbit. Here the final orbit is completely known and the initial orbit is unknown, thus, the 
term "backward” targeting. 

The computer time required to find the optimum transfer is important from a 
practical consideration. It is easy to imagine the large number of computations which 
are performed since the solution is iterative and contains the first and second numerical 
partials of the target parameters with respect to the controls. These numerical partials 
require that the equations of motion (eqs. (16)) be integrated repeatedly. If all six of the 
controls are free, then 12 trajectories are required to compute the first partials and 
61 trajectories for the second partials. Thus, the equations of motion must be integrated 
73 times for a single iteration. In general the required number of trajectories is 
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4n’. 


+ 1) where n is the number of free control parameters. The number of 


2l(n - 2)*. 

iterations necessary to converge the initial guess to the optimum control vector varies 
considerably depending on how good the initial guess is and the sensitivity of the solution. 
However, when the convergence is slow, a large ntunber of trajectories must be inte- 
grated. Although many other calculations are performed, the majority of the machine 
time is spent integrating these trajectories. For this reason the equations of motion 
were expanded in a power series in time and are integrated in this manner. (See the 
appendix.) This approach is much faster than the more standard numerical methods of 
integration. To further decrease the machine time, the trajectories are divided into N 
increments each of which is integrated with a 13 -term power series. Thus, the duration 
of each segment is tj^/N where N = 1,2,. . .,10. As a result of integrating by power 
series, the computer program is not limited by machine time. 


Program Input 

All input to program VTTAP is accomplished by means of a FORTRAN namelist 
DAT. Each of the name list variables is defined in table I. 

The combination of free control variables and desired target parameters is con- 
trolled by an array of 12 integers. NOPT(l) to NOPT(6) correspond to a, /3, 6, 0, 
t| 3 , and Vq. If a 1 is input in NOPT(K), then the Kth control is free to vary. A 0 denotes 
that the Kth control is fixed and not allowed to vary from the Initial guess which is input 
by the GS array. For example, consider the case where the thrust vector is not allowed 
to pitch but remains in a constant inertial direction throughout the bum. Obviously, the 
three angles of,/3,6 (sketch 1) overdefine the problem since only two angles are needed 
to define the inertial direction. To overcome this problem the angle ^ could be held 
constant at 90° allowing or and 6 to function as the right ascension and declination 
of the thrust vector, respectively. Thus, a constant inertial maneuver requires that 
0 = 0 and /3 = 90° throughout the maneuver. The appropriate input for the first six 
integers of NOPT would be NOPT = 1, 0, 1, 0, 1, 1, and GS = 90., 6j, 0., t^^^, 

where the subscript denotes the first guess at the controls. Mathematically, the number 
of control variables is reduced from six to four so that the n of equation (7) is equal 
to four. This is the equivalent of omitting two rows and two columns of the general 
Newton -Raphson matrix of partial derivatives. An interesting point is the number of 
trajectory integrations needed for each iteration. Since only 33 trajectories are needed 
to calculate the partial derivatives instead of the usual 73, this case should use less than 
half the machine time needed for the general case. The constraints imposed on the 
optimum solution are defined by NOPT(7) to NOPT(12). At most VTTAP will consider 
six constraints which would correspond to 1 for NOPT(7) to NOPT(12). However, if 
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NOPT(6 + K) = 0, then the Kth constraint is not invoked. In other words NOPT turns the 
constraints off and on by an input of 0 or 1, respectively. In this manner from one to n 
constraints are considered where n is the number of free controls. In this example 
where /3 and d were fixed the number of constraints must be four or less. The 
specific target parameters are denoted by KOPT, an array of six integers each of which 
have a value from 1 to 5. Table n defines the various options. The values of the target 
parameters chosen are then fixed at the values in array AIN. Therefore, the constraint 
is turned on by NOPT, the target parameter is defined by KOPT, and its value is input 
by AIN. For example, suppose 

NOPT = 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0 

KOPT = 1, 2, 1, 1, 1, 1 

AIN = 20488., 0., 10., 0., 0., 0. 

From NOPT it can be seen that Q!,6,t]j,t'Q are the free controls while |3 and 9 
are fixed. In addition two constraints are considered. KOPT shows that the two target 
parameters are a and i while AIN states that a* = 20488 km and i* = 10°. The 
zeros in AIN have no function since these constraints are turned off by NOPT. However, 
the corresponding input in KOPT is meaningful since the values of the six target para- 
meters chosen from table II are printed as output. Note that one of the target parameters 
in table II is 1/a. This parameter is equivalent to a and provides a smooth transition 
between hyperbolic and elliptical orbits. 

The remaining input parameters of table I need little explanation. The five initial 
orbital elements are defined by the CONI array. The sixth element, true anomaly, is not 
input since it is a control variable. The GL array is similar to GS and contains the ini- 
tial guesses on the Lagrange multipliers. As mentioned previously, the rate of conver- 
gence is highly dependent on the initial guesses at the controls. However, the guesses at 
the multipliers seem to have little effect on the convergence. If GL is not input, VITAP 
fills the array with I's which are acceptable initial values for the Lagrange multipliers. 
The small increments in the control parameters Aof^ used to generate the numerical 
partial derivatives (eqs. (13) and (14)) are defined by the HP array. The maximum allow- 
able step size for the controls during the Newton-Raphson iteration Aof* are defined by 
the VI array. The step size is then computed according to equation (11). Values for both 
HP and VI are built into the program and seldom need changing. However, these values 
may be altered at any time by the appropriate input. The number of segments used to 
integrate the equations of motion by power series is input through NSTEPS. If the total 
burn time t^^ is 1800 seconds and NSTEPS is 6, then VITAP integrates in 300-seconds 
segments with a 13 -term power series, the initial conditions for the second segment being 
the end conditions of the first segment. A reasonable input value for NSTEPS is one 
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which results in about 300- to 400-second segments. The accuracy of the trajectory, 
however, can be improved by increasing NSTEPS at the expense of machine time. The 
three parameters describing the spacecraft are mass, mass -flow rate, and thrust which 
are input according to table I. The gravitational constant of Mars ii is stored in VITAP 
as UMARS = 42828.4 km Vsec^ but can easily be changed by input. Some of the target 
parameters of table II require additional input other than just the value input through the 
AIN array. For example, all the target parameters related to a reference date require 
two Julian dates, the time at which the spacecraft would reach periapsis on the initial 
conic and the time at which the target parameter is to have the given value. These two 
time parameters are PERJD and REFJD, respectively, and need be input only if the time- 
related target parameters are exercised. To constrain the latitude of the landing point, 
the angle 9 must be input through PER. In addition, if the lighting angle G is one of 
the objects of the targeting, then the right ascension and declination of the subsolar point, 
SLONG and SLAT, must also be input. Once all the input is defined, the program VITAP 
uses the initial guesses at the controls and multipliers to start the Newton -Raphson itera- 
tion technique which continues until one of two situations occurs. If the sum of the errors 
(eq. (12)) is less than the value of e which is input through ERR, then the process is 
assumed to have converged at the optimal set of controls and the solution is completed. 

On the other hand, if the number of iterations exceeds the input value of MAXIT, the pro- 
cess is considered nonconvergent and is stopped. 

Finally, the mode of operation must be defined as 1, 2, or 3 by the input parameter 
MODE. The normal targeting mode (MODE = 1) is straight forward. The Initial orbit is 
specified by ao,eo,io>^o>^o along with the desired target parameters. The program 
VITAP finds a control vector 3 which maneuvers the spacecraft from the initial orbit 
to the final orbit which satisfies the target parameters. If there is at least one more 
free control than target parameters, VITAP minimizes the burn time - or equivalently 
the required fuel - by the Newton -Raphson technique (eq. (7)). If the number of target 
parameters is equal to the number of free controls, the solution is a straight forward 
search without minimization (eq. (10)). To target the inclination within bounds, the second 
mode (MODE = 2) is used. This search is similar to mode 1 with the restrictions that the 
inclination of the final orbit must not be specified as a target variable and the number of 
free controls must be greater than the number of target parameters. If the inclination 
constraint option NOPT(9) is not equal to 0, it will be set equal to 0. First the program 
finds the optimum controls which will result in some inclination ij. This inclination is 
then compared to the two boundaries specified by the input parameters BOUND(l) and 
BOUND(2). If the inclination is between these bounds, the solution is complete; if it is 
outside this interval, the program retargets to the closest bound. With the addition of 
another target parameter i, the second solution may no longer be a minimization. 
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The third option (MODE = 3) is more involved than the other two. It has been 
termed backward targeting and finds the best orbit from which to establish a given orbit. 
The only change in input involves the CONI array. Instead of CONI containing the initial 
orbit, for mode 3 it contains the final orbit while the initial orbit is partially defined by 
the target parameter options of table n. Thus, the final orbit is completely defined 
except for the true anomaly and the initial orbit is only partially defined by the target 
parameter constraints. The targeting problem, then, is to find the remaining initial 
orbital elements such that a transfer from this orbit to the given orbit is optimum. It is 
solved with the same mathematics as mode 1 with some changes in procedure. The pro- 
gram internally changes the input and considers the motion from the final orbit to the 
initial orbit. When the motion is reversed the final orbital elements become Bq, 

180*^-io> and ISO^-Wq which is due to the Interchange of the ascending and 

descending nodes. The initial mass of the spacecraft on the initial orbit is known from 
input, but the mass on the final orbit is not known since the optimum bum time has not 
been established. It is obvious that the spacecraft will lose mass during the transfer from 
the initial to the final orbit. Equation (15) can be used to predict the mass in the final 
orbit as mj = mQ + mt^^, where mQ is the initial mass and tjj is a guess at the burn 
time. The initial guesses at the controls are changed internally to a, / 3 , S+et^j, -6, tj^, -Vq 
due to the reversal in motion. The program then uses the corrected orbital elements 
and controls to integrate the equations of motion which results in some initial orbit 
which neither satisfies the constraints nor is optimum. The guesses at the controls 
are then corrected according to the Newton -Raphson matrix, and the iteration procedure 
continued until the solution is obtained. The only difference between this procedure and 
the normal procedure is that the mass mf in the final orbit is corrected for each itera- 
tion according to the change in burn time tj^. Therefore, the solution is complete when 
both the controls and mass have converged. For this reason the backward targeting con- 
verges much slower than the normal targeting. 

Program Output 

The computer outputs for each of the three modes of operation are presented in 
tables in, IV, and V. 

The normal targeting mode (MODE = 1) is presented in table HI. The first output is 
a complete listing of all of the parameters input through the namelist DAT. These param- 
eters have been previously defined in the section entitled "Program Input" and in table I. 
Next is a formal listing of the mode of operation, the initial orbital elements, and a 
description of the free control parameters and the selected target parameters. It can be 
seen that the initial conic is a hyperbola and that the initial guess at the true anomaly is 
-60°. The program will find a subset of four controls ^o!, 6 ,tj 3 ,P'o) such that the final 
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orbit has a prescribed value of a (or l/a), i, and co. Since there are more control 
variables than target parameters, the program will find the controls which minimize 
t]j2. When four controls are varied to find three parameters in the final orbit, the opera- 
tion is said to be a 4 by 3 search. 

The search or iteration is recorded by a selected set of output parameters. The 
first, second, and last iterations are shown in table HI. The control parameters for the 
first iteration which are the initial guesses are shown followed by the initial guesses at 
the Lagrange multipliers. Next, is presented the orbit which resulted from applying the 
initial controls to the initial orbit. The corresponding target parameters are then output 
which are naturally in error since the initial guess at the control variables was not the 
optimum set. The sensed velocity corresponding to t^ is output followed by the error 
in the target parameters. The sensed velocity is the integral of the acceleration due to 
thrust, that is 


AV = 



Thrust 
mQ + mt 


Thrust ^ 

m \ mo / 


(18) 


and the errors in the target parameters are the differences between the desired values 
and the values obtained from the present set of controls. For example, from the input 
it can be seen that the inclination of the final orbit is 35°. This is denoted by the third 
entry in the AIN array. However, the initial controls produced an orbit with an inclina- 
tion of 33.328°. Thus, the third target parameter is in error by 1.672°. Note that three 
of the target parameters have zero error because the second, fifth, and sixth 

target options were not exercised as can be seen from the NOPT array. The errors in 
the target parameters are then used to improve the initial guess at the controls (eq. (7)). 
The corrections to the initial controls are in the next line of printout. Two points should 
be made about these corrections. First, no correction is added to ^ and 9 because 
zeros were input in the NOPT array in the locations corresponding to /3 and 6; thus, 
the program does not allow /3 and 6 to vary but keeps them fixed at the initial values 
which for this iteration were 90° and 0. deg/sec, respectively. The second point is that 
the corrections to the controls were limited by the maximum allowable step size in Uq. 
According to the VI array, the step in true anomaly for one iteration cannot be greater 
than 15°. Since the correction is exactly 15°, this implies that the Newton -Raphson 
technique (eq. (7)) calculated a larger correction and that the corrections were reduced 
according to equation (11). At the end of this block of output is a row labeled 
"Eigenvalues of second partials of augmented function." This row represents the eigen- 
values of the upper left-hand portion of the matrix composed of partial derivatives 


9gi 

(eq. (7)). In other words, they are the eigenvalues of the matrix or 
These values are used as a diagnostic tool. ' 


8^F 

9Qtj 9 Q!j ' 
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The printout of the second iteration follows the same format as the first iteration. 
The control parameters and the Lagrange multipliers have been updated by the correc- 
tions computed in the first iteration and should produce a final orbit that is closer to the 
optimum solution. However, due to the interplay between satisfying both the constraint 
equations and the necessary equations for a minimum, it is not always apparent that the 
succeeding iterations are closer than the preceding ones. For the case being considered 
the second iteration does, however, seem to be closer to the solution. The errors in two 
of the target parameters i,co were reduced while the error in 1/a increased slightly. 
In addition, the sensed velocity was reduced. As before the errors in the target param- 
eters are related to the corrections in the control variables and the iterative process 
continues. 

The Newton-Raphson technique is repeatedly applied until both the constraint equa- 
tions and the necessary equations for a minimum are satisfied (eq. (12)). The 29th 
iteration shows that the process has converged since the errors in the target parameters 
are very small and the computed corrections to the control variables are negligible. The 
optimum set of controls, then, are listed in the first line of printout. They not only pro- 
duce a final orbit with the desired target parameters l/a, i, oj but also establish it with 
a minimum of fuel. Since the eigenvalues of the second partials of augmented function 
are all greater than zero, this is indeed a minimum solution. As a point of interest the 
Newton-Raphson matrix and its inverse (eq. (7)) are output for the last iteration. The 
matrices presented are the full 12 by 12 matrices. The zero rows and columns corre- 
spond to fixed control parameters and to the constraints that were not exercised. 
Accordingly, these matrices are reduced to the proper dimensions during the actual 
computations. The final output is the machine time in seconds required to compute the 
optimum set of control parameters. 

Table IV presents a sample case of targeting such that the inclination is within 
bounds (MODE = 2). The second mode of operation is similar to the normal targeting 
mode and allows for the optimum transfer with the additional constraint that the inclina- 
tion of the final orbit be between an upper and lower bound. The procedure is first to 
disregard the inclination constraint and solve for a set of optimum controls. If these 
controls produce a final orbit that satisfies the inclination inequality, then the solution 
is complete. If this is not the case, the program then determines which inclination 
bound is closest to the present solution and retargets to that value of inclination. This 
procedure is demonstrated in table IV. The transfer is from a hyperbola to an ellipse 
and is a 6 by 4 search. It is necessary that the number of control parameters exceeds 
the number of constraints since the inclination will be added as a constraint if it is not 
within bounds. Also note that inclination is a free variable for the first part of the solu- 
tion and that the final orbit is specified by a, e, co, and Ci. The program proceeds 
to tanget to a, e, co, and Q,. After 12 iterations, the optimum controls are found 
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which produce a final orbit with an inclination of 35.28°. From the input it is seen that 
the constraint on the inclination is 31° = if = 35°. Thus, the inclination constraint was 
not satisfied by targeting to a, e, cu, and O. The program then proceeds to turn 
on the third constraint option (NOPT(9) = 1), define the value of the inclination constraint 
as the closest bound (AIN(3) = 35.), and set the initial guesses at the controls and 
Lagrange multipliers equal to the values obtained in the first search. These changes can 
be seen by comparing the first and second listing of the namelist parameters. The pro- 
gram then proceeds to solve this 6 by 5 search and after 6 iterations converges to the 
optimum set of control parameters. 

The third mode of operation, backward targeting, is presented in table V. This 
mode (MODE = 3) considers the problem of finding the best orbit from which to establish 
a given final orbit. For the particular case shown the final orbit is an ellipse specified 
by the five orbital elements a, e, i, o), and S2. The initial orbit is only partially 
specified by three hyperbolic parameters: the hyperbolic excess velocity, the declination 
of the approach asymptote, and the right ascension of the approach asymptote. These 
three parameters are designated as target parameters. The problem, then, is to deter- 
mine the remaining hyperbolic elements and a set of controls so that when these controls 
are applied to the approach hyperbola the resulting orbit will be the given ellipse. In 
addition the fuel required for the transfer is to be a minimum. The procedure, as out- 
lined previously, is to target backwards from the final ellipse to the initial hyperbola. 

This reversal of motion is handled internally by the program with one exception. The 
initial guess on the true anomaly at the start of the burn Vq is replaced by a guess on 
the true anomaly at the end of the burn. That is, the input value of Uq is 45° which is a 
guess at the true anomaly in the ellipse where the burn will terminate. 

The printout of the first iteration shows the result of reversing the motion and 
applying the intital controls to the ellipse. As would be expected the hyperbolic excess 
velocity, the declination of the asymptote, and the right ascension of the asymptote are all 
in error. These errors are used to correct the control variables. One main difference 
between this mode and the others is the problem associated with the mass of the space- 
craft. In order to integrate the equations of motion from the ellipse to the hyperbola, the 
mass of the spacecraft in the ellipse must be known. The mass, however, is a function of 
the burn time t^, required to establish the ellipse, and this time is not known until the 
solution is complete. Thus, an estimate of the mass in the ellipse is computed as a func- 
tion of the present value of the burn time. This estimate is included in the output under 
the heading "FINAL MASS." Therefore, the iteration process must converge the final 
mass as well as the set of control variables. This additional requirement slows the 
iteration process considerably. For the case considered here the program required 
46 iterations to converge the initial controls to the final set of optimum controls. The 
"best" hyperbola for such a maneuver is given by the third line of printout. 
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CONCLUDING REMARKS 


This report has described a program for determining the minimum -burn -time, 
thrusting, transfer trajectory between two Keplerian orbits. Basically, the method 
described involves the solution of a constrained minimization problem by use of constant 
Lagrange multipliers and the Newton-Raphson iteration technique. The model of the 
transfer trajectory allows for the thrust vector to rotate in a plane with a constant pitch 
rate. In all, six control parameters characterize the transfer trajectory; thus, as many 
as six constraints can be imposed on the final orbit. If the number of constraints is less 
than six, the program solves for the set of controls which minimizes the burn time or 
equivalently the required fuel. The option is available to fix any of the six control vari- 
ables and allows various guidance laws to be investigated such as a constant inertial burn. 
This option together with the twenty different constraints from which to choose allows a 
considerable amount of flexibility. In addition the program VITAP operates in three 
modes. The first mode is the normal optimum transfer, the second mode is similar to 
the first with the additional constraint that the inclination of the final orbit be between an 
upper and lower bound, and the third mode targets backwards. This mode considers the 
problem of finding the best hyperbola from which to establish a given ellipse. These 
three modes of operation have been exercised extensively for representive 1975 Mars 
orbital insertion maneuvers. 

The program uses only necessary conditions for a minimum. Accordingly, there is 
no guarantee that the solution obtained will be a minimum. If, however, the eigenvalues 
of the matrix corresponding to the second partials of the augmented fmction are all posi- 
tive, then the solution is indeed a local minimum. This particular problem of using only 
necessary conditions offers no great difficulty when there is a judicious choice of trajec- 
tory constraints and reasonable initial guesses on the control parameters. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 5, 1971. 
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APPENDIX 


POWER SERIES EXPANSION OF EQUATIONS OF MOTION 


The most time consiiming computation in the Viking Targeting Analysis Program 
is the numerical integration of the equations of motion. Since numerical partial deriva- 
tives are used in the Newton-Raphson technique, the equations of motion are integrated 
repeatedly. For this reason the equations of motion are numerically integrated by power 
series which are considerably faster than the more standard numerical methods. 

The power series solution of the two body equations of motion have been solved by 
Schanzle (ref. 11). The governing differential equations for two body motion are 


where 


r3 




(Ala) 

(Alb) 

(Ale) 


(A2) 


and iJ. is the gravitational constant. The trajectory model considered for the Mars 
maneuvers requires the addition of acceleration terms which are due to the low thrust 
(nonimpulsive) engine. Therefore, the equations of motion for the finite burn are 
(eqs. (16)) 


X = — ^ + a(t) [cos (6+0t) cos a - sin (6+0t) sin a cos 
r3 

(A3a) 

y = — ^ + a(t) [cos (6+0t) sin a + sin (6+0t) cos or sin /£] 

r<3 

(A3b) 

z = + a(t) [sin (6+0t) sin 

r3 

(A3c) 

where a(t) is the magnitude of the acceleration and is given by 


a(t)= Thrust 
mo + mt 

(A4) 
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APPENDIX - Continued 


The following definitions will simplify the solution of the equations of motion: 



Ax(t) = a(t) [cos (6+0t) cos a - sin (6+^t) sin a cos ^ 
Ay(t) = a(t)[cos (6+0t) sin a + sin (6+0t) cos a cos ^ 
A 2 (t) = a(t) [sin (6+0t) sin ^ 

Therefore, equations (A3) can be written as 
X = -jtibx + Ax 
y = -jLiby + Ay 
z = -jubz + A 2 

for which the assumed solution is 

i=oo 

X = ^ Xjt^ 

i=0 

i=oo 

y = I yi‘* 

i=0 

i=oo 

i=0 

i=oo 

y = I yi‘‘ 

i=0 

i=oo 

i=0 

i=oo 

■^x = ^ •^Xji*^ 

i=0 


(A5) 

(A6a) 

(A6b) 

(A6c) 

(A7a) 

(A7b) 

(A7c) 

(A8a) 

(A8b) 

(A8c) 

(A8d) 

(A8e) 

(A8f) 
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APPENDIX — Continued 


i=co 



(A8g) 

i=0 


i=oo 


1 

(A8h) 

i=0 



The three position coordinates x,y,z along with the radius r, the parameter b, 
and the acceleration components Ax,Ay,Az are assumed to be represented in some 
neighborhood of assumed values by Taylor series expansions. Taylor series or power 
series will be used interchangeably in this appendix which can be justified by Theorem 39, 
page 354 of reference 9. The sufficient conditions for convergence of the particular 
series are covered, at least for the two body case, by Schanzle (ref. 11) and should be 
easily extendable to the present case of a constant thrust and constant mass -flow -rate 
burn. In fact, the convergence of the series for the acceleration terms of equations (A3a) 
to (A3c) is trivial, being the multiplication of the series for the sine and cosine of ^t, 
which is absolutely convergent everjrwhere, by the series for a(t), which is absolutely 
convergent for t < mo/m. (See ref. 9.) Of course, t cannot get as large as mQ/m 
which corresponds to the spacecraft mass becoming zero. 

In order to obtain the solution to equations (A7a) to (A7c) it is necessary to evaluate 
the coefficients of equations (A8a) to (A8h). The three series for the acceleration terms 
(eqs. (A8f) to (A8h)) are determined from equations (A6) while the series for x, y, z, 
r, and b are determined by means of recursive equations. The following general for- 
mula will aid in this development (ref. 12), 



i=oo/k=i 

E I Pk9l-k)‘‘ 

i=0\k=0 


(A9) 


where p and q are the coefficients of the general power series. The recursive equa- 
tion for X can be developed as follows. Differentiating equation (A8a) twice gives 

i=oo 

X = ^ i(i-l)xjt^‘^ 

and then substituting into equation (A7a) along with equations (A8e) and (A8f) yields 
i=oo /i=oo \A~°° \ 

i=0 \i=0 /\i=0 / i=0 
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APPENDIX — Continued 


Changing indices on the first summation gives 


1 = 00 


^l=co \ / i=co 


1=00 


^ (i+l)(i+2)xi+2t^ = -Ml ^ hni J XitM + ^ Ax,it^ 


i=0 


\i=0 /\i=0 


i=0 


(AlO) 


or by use of equation (A9) it may be rewritten as 

i=oo i=oo k=i i=oo 

^ (i+l)(i+2)xi+2t^ = 2 ^ Ax,it^ 


(All) 


i=0 


i=0 k=0 


i=0 


Equating the coefficients of the nth power of t gives the recursive equation for the coef- 
ficients of the X series, that is 


k=n 

(n+l)(n+2)xn+2 = "M ^ + A^^n 

k=0 


(A12) 


The recursive equations for y and z are similar. Therefore, the three recursive 
equations for the coefficients of the x, y, and z series are 


Xn+2 


1 

(n+l)(n+2) 



k=n 

^ Vn-k + 
k=0 



yn+2 


1 

(n+l)(n+2) 



k=n 


X Vn-k 

k=0 


+ 



(A13a) 


(A13b) 


^n+2 


1 

(n+l)(n+2) 



+ 



(A13c) 


where Xq, yQ, zq are the components of its initial position and Xj, yj, z^, are the 
components of the initial velocity. 

The recursive equation for r is found by differentiating equation (A2), that is 


ri* = XX + yy + zz 

and by substituting equation (A8a) to (A8d) to yield 
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APPENDIX — Continued 


1=00 \ /i=oo \ /i=°o > 


i=0 / \i=0 


/ \i=0 


1=00 


^ (i+l)xi+it^ 
.i=0 


/i=oo 

^ I v.‘‘l 

\l=0 


1=00 




i=0 


' 1=00 


\l 


I 


a=0 /Li=0 


i=oo 

^ (i+l)zi+it^ 


Multiplying the series together by equation (A9) and equating the nth power of t gives 


k=n 

X (^+l)rk+l^n-k 

k=0 


k=n 

X (k+l)(xk+lXn_k + Yk+lYn-k + ^k+l^n-k) 
k=0 


Removing the nth term from the summations produces the recursive equation for the coef- 
ficient of the r series, that is 

"‘"+1 " ^n+l^o) 

k=n-l 

+ ^ (k+l)(xk+iXn_k + yk+lVn-k + ^k+l^n-k " ^k+l^n-k)] (A14) 

k=0 ^ 


where 

^o= (xo + yi + z?)^^^ 

(Vi + yoyi + zqZi) 

The recursive equation for b is found by differentiating equation (A5), which gives 
rb = -3br 


In a manner similar to that used to develop the recursive equation for r one obtains the 
recursive equation for b as 


^n+1 “ 


1 

(n+l)ro 


k=n-l 


-3(n+l)rn+lbo - ^ (^^+l)(3rk+lt'n-k + Vl^n-k) 


k=0 


(A15) 
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APPENDIX - Continued 


where 


bo = 


and 


bi = 


-Sborj^ 


Through equations (A13), (A14), and (A15), the coefficients of the assumed series solution 
to the differential equations of motion can be found which lead directly to the state of the 
spacecraft at time t by evaluation of the power series. The only remaining development 
is the series expression for the three acceleration terms. 


The power series expansions of the three acceleration functions are most easily 
accomplished by considering the product of two series, the first of which is a(t). Expan- 
sion of the magnitude of the acceleration (eq. (A4)) yields 


a(t) = 


Thrust 
mo + iht 


= Thrust 



(A16) 


The second series are the direction cosines which can be rewritten by use of the general 
trigonometry summation formulas as 


[cos (&f0t) cos a - sin (6+0t) sin a cos /^ = Fi cos 0t + F2 sin 6 t 
[cos (6f0t) sin a + sin (6+0t) cos O! cos |3| = F3 cos 6t + F4 sin 6t 
[sin (6+0t) sin ^ = F5 cos 0t + F0 sin 0t 


where Fj^ is a constant and is given by 


Fj^ = cos a cos 6 - cos /3 sin a sin 6 
F2 = -cos a sin 6 - cos |3 sin a cos 6 
F3 = sin or cos 6 + cos /3 cos a sin 6 
F4 = -sin a sin 6 + cos ^ cos a cos 6 
Fg = sin j3 sin 6 
Fg = sin /3 cos 6 

The definitions for sin 6 t and cos 0t can be given in terms of a power series as 
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APPENDIX — Concluded 


• M 05 t 5 

sin 0t = 0t — — + 

(J • 


I — V. — 

i=0 


^ 1 ©2^2 04t4 

COS et = 1 — — + 


l=c» 


4 ^ cos(^^)eHi 

Z — ly— 


i=0 


Therefore the three acceleration functions can be written as 

-1 r 


A^(t) = 


Av(t) = 


1=00 


Thrust 




i=oo 

i =0 “^o J [i=o L 


i+1 




il 


+ F. 


i: 


1=00 


Thrust 


(-l)ifh¥ 


i=0 “o 


i+1 


1=00 


i=0 


- ^ -(f 

Fs — n — +F 4 — - — 


AzW - 


i=oo . • • 1 1=°® 

Thrust y H)^ J y 

i4 <*' liio 


Sln(f>‘ cos(f)«f 


Fs tf^+Fe „ 


1: 


Multiplication of the two series by the general formula of equation (A9) yields the desired 
power series expansion for the three acceleration functions as follows: 


i=oo k=i 


Ax(t) = Thrust M 2 


i=0 k=OL 


Fi sin(^y^ + Fg cos(|^)0: 


kl 


i=oo |k=i 

Ay(t) = Thrust ^ \ ^ 
i=0 |^k=0 

i=oo |k=i 

AgCt) = Thrust ^ \ ^ 

i=0 I k=0 


F3 sin(^)^ + F4 cos(|^) 0 : 


kl 


F5 sin(|^)0k + Fg cos(^)e^ 
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TABLE I 


DEFINITION OF INPUT PARAMETERS FOR PROGRAM VITAP 


Program 

symbol 

Mathematical 

symbol 

Dimensions 

Units 

Definition 

NOPT 


12 

None 

Integer array denoting free control variables and selected target 
parameters, N0PT(1 to 6 ) corresponds to 
NOPT® = 0 for Kth control fixed and NOPT® = 1 for Kth 
control free where K = 1,2,. . ., 6 ; N0PT(7 to 12) corresponds 
to • .,^ 0 ; N 0 PT(K 46 ) = 0 for Kth constraint omitted; 

NOPT(K+ 6 ) = 1 for the Kth constraint considered 

KOPT 


6 

1 

None 

Integer array denoting the specific target parameter chosen (see 
table II) 

AIN 


6 

1 

km, sec, 
or deg 

Array of values for target parameters denoted by KOPT 

CONI 

j ^ 0 »® 0 ’ ^ 0 » ^ 0 ? ^0 

! 5 

km or deg 

j Orbital elements of initial orbit; ag < 0 for hyperbola 

GS 

1 

6 

deg or sec 

Initial values (guesses) of controls of,i3,6,0,t^,i/Q (see 
sketch 1 ); these values will vary or remain fixed depending on 
NOPTd to 6 ) 

GL 

: M 

6 

j 

Initial guesses on Lagrange multipliers; if not input, 
GL(1 to 6 ) = 1. 

HP 

Ao?i 

6 

i deg or sec 

Increments of controls for numerical partial derivatives; if not 
input, HP(1 to 6 ) = 0.6, 0.6, 0.6, 0.006, 1, 0.6. 

VI 

Aq?* 

1 

6 

deg or sec 

’ Maximum allowable step size for controls during Newton -Raphson 

iteration; if not input, Vl(l to 6 ) = 15., 15., 15., 0.0015, 50., 15. 

NSTEPS 

1 ; 

1 

None 

Integer denoting number of segments used for power series 
solution to equations of motion; if not input, NSTEPS = 10. 

MASS 

• 1^0 

1 _ . . . 

1 

, kg 

1 

Initial mass of spacecraft 



TABLE I - Concluded 

DEFINITION OF INPUT PARAMETERS FOR PROGRAM VITAP 


Program 

symbol 

Mathematical 

symbol 

Dimensions 

Units 

Definition 

DMASS 

m 

1 

kg/sec 

Mass -flow rate 

THR 


1 

kN 

Thrust of the spacecraft propulsion system 

PERJD 


1 

days 

Julian date of periapsis passage on initial conic 

REFJD 


1 

days 

Reference Julian date for constraints of longitude, latitude, and 
true anomaly at a reference time; PERJD and REFJD need not be 
input if these constraints are not used 

PER 

e 

1 

deg 

Angular distance from periapsis to landing point (see sketch b) 

SLAT 


1 

deg 

Declination of subsolar point in areocentric equatorial coordinate 
system 

SLON 


1 

deg 

Right ascension of subsolar point; Input SLAT and SLONG only if 
KOPT(5) = 5 

ERR 

e 

1 

None 

Newton -Raphson convergence criteria (see eq. (12)); if not input, 
ERR = 10"^ 

MAXIT 


1 

None 

Integer denoting maximum number of iterations allowed; if not 
input, MAXIT = 50. 

UMARS 1 


1 

kmVsec^ 

Mars gravitational constant; if not input, UMARS = 42828.4 

MODE 


1 

None 

Program mode: 

1 - normal forward targeting; 

2 - forward targeting, inclination within bounds; 

3 - backward targeting 

BOUND 


2 

deg 

Bounds on inclination for MODE = 2; BOUND (1) is lower bound 



TABLE n 


TARGET PARAMETER REQUEST KEYS 


Input 

value 

Input parameter 

KOPT(l) 

K0PT(2) 

K0PT(3) 

K0PT(4) 

K0PT(5) 

K0PT(6) 

1 

a 

e 

i 

0) 

Cl 

f 

! 

2 

1 

a 

ra 

^-Latitude at 
reference date 

^Longitude at 
reference date 

^Longitude at 
reference date 

^True anomaly at ' 
reference date 

3 



B-t 

(see sketch 3) 

^Latitude at 
reference date 

^^Latitude at 
reference date 

^Longitude at 
reference date 

4 

1 

Orbital period 

1 

1 

1 

B-R 

(see sketch 3) 


t>Declination 
of S 

'^Right ascension 
of S 


i 

i 

1 5 

! 

i 

P 

’ V 

I ’ 00 



^Latitude of 
landing point, 0 

^Lighting angle at 
landing point, G 



^ Value of PERJD and REFJD must be input (see table I). 

^S H Incoming hyperbolic asymptote (see sketch 3). 
c Value of PER must be input (see table I), 

^ Value of PER, SLAT, and SLON must be input (see table I). 



TABLE III 


SAMPLE OUTPUT FOR NORMAL TARGETING MODE 
[Mode = l] 

$DAT 

NOPT * I, C, 1, 0, 1, 1, U 0, 1, I, C, C, 

KOPT = 2, 1, 1, 1, 1, I, 

AIN = D.4€e878024932E-C^, O.Ct 0,35E+02» 0.13C38E+03, 0*0t 0,0» 

CONI = -0.AC39C96E+C4, C.197CAE+01, 0.3531E+C2* 0. 1290692335 1 lE+03t 

0.3C3C7I8220A6E+C3, 

GS = -0.5E4C2» C,9E + 02» 0.2E + 02t O.C» 0.24E+0A, -C.6E-I-02, 

GL = 0.1E4C1» C.lE+01, O.IE+Cl, 0,1E+01, O.IE+Cl, C.lE^Ol, 

HP » 0.6E + CCf 0.6E+0C , 0.6E+00, C.6E-03f O.IE + Ol t 0 *68 + 001 

VI = 0*15E402, 0.15E+C2t C.15E+02, C.15E-02, C.5E+02, 0.15E+02, 

NSTEPS = 8, 

MASS = 0.320695E+CA, 

DMASS = -O.A172E+OC, 

THR = 0.133A5E4C1, 

PERJD = If 

REFJO * 0*0i 

PER =1, 

SLAT = 1, 

SLON * 1, 

ERR = 0.1E-C7, 

MAXIT = 50, 

UMARS = 0.^282848+05, 

MODE = I, 

BOUND = 0.0, 0. 188+03, 

SEND 


CO 



TABLE ni - Continued 


VIKING TARGETING ANALYSIS PROGRAM IVITAP) 


♦♦♦♦♦NORMAL TARGETING (MOOE=l)^^^^^ 
INITIAL CONIC 


SMA 

-4939.396C ECC 1.S7C4000 

INC 

35.310000 

PER 129.06923 

NOD 303.07182 

TAN -60.000000 

THE 

CONTROL VAPIAfiLES ARE ALPHA 

DELTA 

TBURN 

TASTART 



THE 

TARGET VARIABLES ARE 1/SMA 

INC 

PERI 





ITERATION 1 






CONTROL PARAMETERS 
ALPHA 

-50.C00CGCC000 

BETA 

90.CCOQCCCOOO 

DELTA 

20.CCCCC0C0CC 

THDOT 

0. 

TBURN 

24OO.0OCC00OO 

TASTART 

-60.0000000000 

LAGRANGE MULTIPLIERS 
l.OOOOOOCCOOO 

0. 

l.CCCCCCCCCCC 

1.00000000000 

C. 

0. 

ORBIT OBTAINED FROM PRESENT CCNTPCLS 
SMA 68636.103 ECC .93C9C775 

INC 33.327990 

PER 139,78924 

NOD 307.84094 

TAN 37.890110 

TARGET PARAMETERS 
1/SMA 

1.456959179706E-05 

ECC 

.930907746419 

INC 

33.3279901131 

PERI 

139.789236539 

NODE 

307.840943497 

TAN 

37.8901097288 

DELTA V 1.23!5143 






ERRORS IN TARGET VARIABLES 

1/SMA ECC 

3.431821C69614E-05 0. 

INC 

1.672CC986694 

PERI 

-9.40923653906 

NODE 

0. 

TAN 

0. 

CORRECTIONS TO CCNTRGL VARIABLES 

ALPHA BETA 

-9.91246578435 0. 

DELTA 

6.1C646498221 

THDOT 

0. 

TBURN 

-17.66C69281C5 

TASTART 

-15.0000000000 


EIGENVALUES OF SECOND PARTIALS OF AUGMENTED FLNCTICN 

7.94149282?7D+01 2. 00004378900*00 -2. 98 1 5 233394D+C 1 -5. 51 829307240*01 



TABLE III - Continued 


mRATIC^ 2 


CONTROL PARAMETERS 
ALPHA 

-59.9124657843 


BETA 

90.CCOCCOOOOO 

delta 

26. 1064649822 

LAGRANGE MULTIPLIERS 
-564455.121282 


0. 

8.40845983881 

ORBIT OBTAINED FRCM 
SMA 132642.89 

PRESENT CONTROLS 
ECC .96780569 

INC 36.227467 

TARGET PARAMETERS 
1/SMA 

7.539C39663114E 

-C6 

ECC 

.967805693215 

INC 

36.2274669615 


THDOT 

TBURN 

2382.33930719 

TASTART 

-75.000C000000 

.7661C8102041 

C. 

0. 

138.67499 

NOD 316.23571 

TAN 351.64946 

PERI 

138.674993190 

NODE 

316.2357C5212 

TAN 

351.649461687 


DELTA V 1.22^1C60 


ERRORS IN TARGET VARIABLES 
1/SMA 

A.13A876263009E-C5 0. 


ECC 


INC 

-l.227A669tlA8 


PERI 

-8.29A99319019 


NODE 


C. 


TAN 


CORRECTICNS TO CONTROL VARIABLES 

ALPHA BETA 

2.E572233C499 0, 


DELTA 

6.77719297AI89E-C2 0, 


THDOT 


TBUPN 

-50*aOOCQCOCQQ 


TASTART 

•1A261C85E923 


EIGENVALLES OF SECOND PARTIALS OF AUGMENTED FUNCTICN 

1.26 52278 0 2ED + 03 1.95 3534 185 6D«-02 1.9 999762086C+C0 -1. 1 3440 34257D+02 


ITERATION 29 


CONTROL PARAMETERS 
ALPHA 

-25.6732430975 

BETA 

9C.COCCCOCOOO 

delta 

22.5386393285 


THDOT 

0. 


TBURN 

2453.84C69840 

TASTART 

-63.8506518549 

LAGRANGE MULTIPLIERS 
-47633505555.7 

0. 

5595C. 9037335 


101.077861597 


C. 

0. 

ORBIT OBTAINED FPCH PRESENT CCNTRClS 
SMA 20455. OOC ECC .76436452 

INC 35.000000 

PER 

130.38000 

NOD 

303.97940 TAN 

40.321847 

TARGET PARAMETERS 
1/SMA 

4.888780249287E-C5 

ECC 

.764364521465 

INC 

35.COOOOOCCCO 


PERI 

130.380000000 


NODE 

303.9794C0569 

TAN 

40.3218474563 

DELTA V 1.27C5839 








ERRORS IN TARGET VARIABLES 

1/SMA ECC 

3.267785347871E-16 0. 

INC 

-2. 7284041O5319E- 

12 

PERI 

1 .909938873723E' 

-11 

NODE 

0. 

TAN 

0. 

CORRECTICNS TO CONTROL VARIABLES 

ALPHA BETA 

3.012653474970E-11 0. 

DELTA 

6.C69949388672E- 

11 

THDOT 

0. 


TBURN 

3.215712718649E-09 « 

TASTART 

-2.175990898215E-10 


EIGENVALLES OF SECOND PARTIALS OF AUGMENTED FUNCTICN 
CO 2.A3A 11942530+07 1.28 835964780+07 3.83341670410+06 1. 96 A85 702520+00 

CO 


o 


TABLE III - Concluded 


NEWTON R4PHSCN MATRIX 
1.267E+C7 0. 

5.835E+C5 

0, 

-2.079E+C3 

-7.553E406 

4.504E-06 

0. 

3.904E400 -2.946E401 

0« 

-2.678E-C1 

0. 

0. 

c. 

0, 

0 • 

0. 

0 « 0 « 

5.835E+C5 

0. 

l.AllE+07 

0. 

-1.9C5E4C3 

-6.888E+C6 

-e.788E-06 

0. 

-7.436E40C -4.256E401 

0. 

0. 

0. 

-2.878E+05 

C. 

0. 

0. 

0. 

0. c . 

-2.079E+03 

0. 

-1.9C5E+C3 

C* 

3,29CE4CC 

4,330E+03 

U020E-07 

0. 

-4.834E-04 -5.321E-03 

-7.553E+C6 

3. 

-6.e88E4C6 

C. 

4.330E+03 

1.A28E+07 

-1.406E-O6 

0. 

-1.225E400 -2.549E401 

4.504E-C6 

0. 

-8.78fE-C6 

0. 

1.02CE-C7 

-1.486E-C6 

0. 

0. 

0. C. 

0. 

3. 

0. 

0. 

C. 

C. 

0. 

0. 

0 . c. 

3.9C4E+C0 

0. 

-7.436E+C0 

0. 

-4.83AE-04 

-1.225E400 

0. 

0. 

0. c. 

-2.946E+C1 

3. 

-A.256E401 

0. 

-5.221E-03 

-2.549E+01 

0. 

0. 

0 . c . 

0. 

0. 

0. 

0. 

C. 

0. 

0* 

0. 

0 . c. 

0. 

0. 

0* 

c. 

C • 

0. 

0. 

0* 

0. c. 

INVERSE 

4.539D-C9 

0* 

4.479D-C9 

0. 

-U93C0-16 

-1.2730-08 

4.4290402 

0. 

1.0440-Cl -1.4730-02 

0. 

0. 

0. 

0* 

0. 

C. 

0. 

0. 

0. 0. 

4.479D-C9 

0. 

4.42CD-C9 

0. 

-1.9C50-16 

-1.256C-C8 

2.5380401 

0. 

-8.1950-C2 -5.6220-03 

0. 

0. 

0. 

0* 

C. 

c. 

0. 

0. 

0. 0. 

-1.930D-16 

0* 

-1.905D-16 

0. 

8.21CD-24 

5. 4130-16 

9.7470406 

0. 

-1.1400401 -2.0600-02 

-U273D-C8 

0. 

-1.256C-C8 

0. 

5.4130-16 

3. 5690-08 

-2.5890403 

0. 

1.850C-C2 -1. 2810-02 

4.^290*C2 

0 . 

2.53804CI 

c. 

9.7470+C6 

-2.5890403 

-1.9160414 

0. 

-3.2770408 9.4190407 

0* 

0 . 

0. 

0« 

0. 

0. 

0. 

0. 

C. C. 

1.04^D-C1 

0, 

-8.195C-02 

0. 

-UU0D4CI 

1.8500-02 

-3.2770408 

0. 

-2.195C4C5 1.0190404 

-1.4730-C2 

0. 

-5.622C-03 

0. 

-2.C6CC-02 

-1.2810-02 

9.4180407 

0. 

1.C19D4C4 -1.7920403 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 

0. 

3. 

0. 

0. 

0. 

c. 

0. 

0. 

0 . c. 


TIME FOR THIS C^SE 


22.7A0 


oooooooooooo OOOOOOOOOOOO 



TABLE IV 

SAMPLE OUTPUT FOR TARGETING WITHIN BOUNDS 

[Mode = 2] 

iDAT 

NOPT = 1, 1, I, 1, 1, 1, 1» I, 0, It It C, 

KOPT = 1, 1, 1, It 1, 1, 

A!N = 0.20455E+C5, C.776t 2S3723t5E*CC, C^C, C. nC36E<-03 1 0. 30327E+03t 

0 . 0 , 

CONI = -0.4S39096E+CA, 0.19704E+01, C.3531E+C2, 0.1290692335116+03, 

0.303071822048E+03, 

GS = -0.55E+C2, C.AE+C2, C.5E+01, 0.25E-C1, C.2^E+C^t -0.6E+02, 

GL = O.IE+Cl, O.IE+Olt O.IE+Ol, C.IE+Ol, O.IE+Cl, O.IE+Cl, 

HP = 0.6E4CC, C.6E+CC, C.6E+C0, C.6E-C3, O.lE^Cl, C.6E+00, 

VI = 0.15E+02, 0.15E+C2, C.15E+02, C.15E-02, C.5E+02, 0.15E+C2, 

NSTEPS = 8, 

MASS = 0.32C695E+CA, 

OMASS = -0.4772E+C0, 

THR = 0.13345E+C1, 

PERJD = It 

REFJO ^ 0.0, 

PER =1, 

SLAT = I, 

SLON = I, 

ERR = 0.1E-C7, 

MAXIT = 50, 

UMARS = 0.42E284E+C5, 

MODE = 2, 

BOUND = 0.31E+C2, C.35E+02, 


SEND 


TABLE IV - Continued 


VIKING TARGETING ANALYSIS PROGRAM (VITAP) 
♦♦♦♦♦INCLIKAT ICN BETWEEN ECUNCS ( MC C E= 2 M ♦ 


INITIAL CCMC 
SMA -AS39,096C ECC 

1.97CACCC 

INC 

35.21CCCO 

PEP 

129.06923 

NCO 303.07182 

TAN -60.003000 

THE CCNTPOl VflPIABLES ARE 

ALPHA 

EETA 

DELTA 

THCCT 

TBURN 

TASTART 


THE TARGET VARIABLES ARE 

SMA 

ECC 

PERI 

NCOE 





ITERATICN 1 






CONTROL PARAMETERS 
alpha 

-55.CCOOCCCOOO 

EETA 

AC.CCCCCOOOOO 

DELTA 

5.CCOOCCOCCCC 

ThOOT 

2.5COOCOOOOOOOE-02 

TBURN 

2AOO.OOCCOCCO 

TASTART 

-6O.00000C00OC 

LAGRANGE MULTIPLIERS 
l.CCCOCCCCOOO 

l.CCCCGCCCOOO 

C. 

1.00000000000 

l.CCCOCCCOCCO 

0. 

CRBIT OBTAINED FPCM 
SMA 193A9.912 

PRESENT CCNTROLS 
ECC .76A66673 INC 

3A.806A52 

PER 128.75868 NOD 

303.52583 

TAN 50.136035 

TARGET PARAMETERS 
SMA 

193A9.912CC75 

ECC 

.7^ A666731799 

INC 

3A.806A515029 

PERI 

128.758679590 

NODE 

303.525826C65 

TAN 

50. 13603A902A 

DELTA V 1.23'^IA3 






ERRORS IN TARGET VAR 
SMA 

11C5. 03799253 

lABLES 

ECC 

1.19626A056613E-02 

INC 

0. 

PERI 

1.621320A1005 

NODE 

-.2558266CA55A 

TAN 

0. 

CORRECTICNS TO CCNTROL VARIABLES 

ALPHA BETA 

-11.8 7A68 6C926 - 1 5 . GCCO QC CCOO 

DELTA 

I1.A010806779 

TFDOT 

-3.259386365AC8E-0A 

TBURN 

-15.2568751A92 

TASTART 

-2.43889172648 

EIGENVALUES OF SECOND PARTIALS OF AUGMENTED 
2.3 1589A77 3 8D+11 8. A 1 282603C ID* OA 

FUNCTION 

5.62A59016A8C+0A 

2.02037A6052D+0A 

1.766CA10A880+03 

2. 1777135A99D+C 



TABLE IV “ Continued 


ITERATICN 2 


CONTROL PARAMETERS 

ALPHA 

BETA 


DELTA 

THOOT 

TBURN 

TASTART 

-66.874666C926 

25.CCCCCCCOOO 


16.401C8C6779 

2 .467406136346E-02 

2384.74312465 

-62.4388917265 

LAGRANGE MLLTIfLIERS 

3.11933539515 

-10592. 155C964 


C . 

42. 1739821634 

-9.18282948849 

0. 

ORBIT OBTAINED FPCM PRESENT CCNTRCLS 
SMA 20204.416 ECC .71402282 

INC 

36.23^337 

PER 131.02962 NOO 

302.30496 

TAN 43.941272 

TARGET PARAMETERS 

SMA 

ECC 


INC 

PERI 

NODE 

TAN 

20204.4162659 

.774022811646 


36.2393371766 

131.029617172 

302.304960945 

43.9412715869 

DELTA V 1.2256561 

ERRORS IN TARGET VARIABLES 

SMA ECC 


INC 

PERI 

NODE 

TAN 

25C.5837241C8 

2.6C6554719133E- 

C3 

C. 

-.649617171922 

.965039054750 

0. 

CORRECTICNS TO CCNTRCL 
ALPHA 

VARIABLES 

BETA 


DELTA 

THOOT 

TBURN 

TASTART 

15.000000COOO 

7.4C529249771 


-13.239C6C7551 

-7.521978869204E-04 

-10.2106029132 

.669162760881 

EIGENVALLES OF SECOND 
7.5295662985D+11 

PARTIALS OF AUGMENTED 
3.09:11345920+05 

FUNCTICN 

2.1743498928D+C5 

9.3584733411D+04 

7.8486027692D+03 

2. 61596658610+00 

ITERATION 12 
CONTROL PARAMETERS 

ALPHA 

BETA 


DELTA 

THOOT 

TBURN 

TASTART 

-58.3595295367 

34.947603C320 


8.89468323431 

2.40746244642CE-02 

2371.64763119 

-61.7465949681 

LAGRANGE MULTIPLIERS 

94.2815423607 

517787.364574 


C. 

182.143038384 

-3404.67367318 

0. 


ORBIT OBTAINED FPCH PRESENT CCNTPOLS 

SMA 20^55,00C ECC .77662S37 INC 35.281270 PER 130.38000 NOD 303.27000 TAN 44.620900 

TARGET PARAMETERS 

SMA ECC INC PERI NODE TAN 

2O455.OQC00OO .716629372365 35.2812702454 130.380000000 303. 2700C0CC0 44.6208995264 

DELTA V 1.2172220 

ERRORS IN TARGET VARIABLES 

SMA ECC INC PERI NODE TAN 

-4.656612313077E-10 -3. 5 f 2713678801E- 15 0. -9.094947017729E-13 C. 0. 

CORRECTICNS TO CCNTRCL VARIABLES 

ALPHA BETA DELTA THOOT TBURN TASTART 

-7. 3813986 30253E-C9 -6. 76 50 1 1 56 E 7C 8E-C9 6. 09 248 1 265466E-09 7. 8852682365 81 E- 14 9. 678 59936565 2E- 12 -1.63799593386 lE-10 

EIGENVALLES OF SECCND PARTIALS OF AUGMENTED FUNCTION 

2. 50717982620+13 9. 759 666 5546D+06 6 .965C75902 lC+06 2. 5602547370D+06 2.4947361034D+C5 2. 23577477660+Cl 


CO 


*♦♦♦♦ INCLINATKN NOT WITHIN BCUNOS - BEGIN TARGETING TC NEAREST BOUND***** 





TABLE IV “ Continued 

$DAT 

NOPT = 1, 1, It If It If If If If If If C, 

KOPT = If 1, I, If If It 

AIN = 0.2CA55E + C5f 0. 7 -J 66 2<5 3 7 2 36 5E +CC , 0.35E + 02f 0.13038E + 03» 

0.3C227E+C3, 0.0* 

CCNI = -tO. 4C39C96E + 04, 0.1970AE + 01, 0.3531E+02f 0 . 1 2 90692335 1 lE + 03 , 

3*3C3C71822CAeE+C3f 

GS = -0.5 63595 295 A40AlE+02» 0 . 3A94 760 3025 2C9E +02 , C . 8 89^683 2A09062E+0 1 , 

3.2AC7A62AA6A28E-01, 0 . 2 37 1 6A 763 119 1 lE+C A » -C . 6 1 7A6594 968311 E+02 , 

GL = 0.g<2815A2ACC895E+C2, C . 5 1 77 6 736 129 5 9 2 E + C 6 , C.Ct 

0. 1821^3038516596+63, - 0. 3 A046 736 587 5 C 3 E+C A , C.O, 

HP = 0.6E+00, 0 *66+00, 0.6E+00, C.6E-03, O.lE+01, C*6E+C0f 

VI = 0.15E + C2, C.15E + C2, 0.156+02, 0.15E-02, C.5E+02, 0.15E+-C2, 

NSTEPS = 8, 

HASS = 0.32C695E+CA, 

DMASS = -0.A772E+00, 

THR = 0.133A5E+C1, 

PERJD = I, 

REFJO = 0.0, 

PER =1, 

SLAT = I, 

SLON = I, 

ERR = 0.1E-C7, 

MAXIT = 50, 

UMARS = 0.A2828AE+C5, 

MODE = If 

BOUND = 0. 3 1E+C2, C.35E+02, 


^ ■ ' 


SEND 



TABLE IV - Continued 




VIKING TARGETING 

ANALYSIS 

PROGRAM 

(VITAP) 


♦♦♦♦♦NOFFAL T4PCETING (M0DE= 

1 )♦♦♦♦♦ 






IMITIAL CCMC 







SMA -A939.0960 ECC 1. 

S7C40CC 

INC 35.310000 

PEP 129. 

06923 

NOD 303.07182 

TAN -61.7A65S5 

THE CCNTFCL VARIABLES ARE 

ALPHA 

BETA DELTA 

THOCT 

TBURN 

TASTART 


THE TARGET VARIABLES ARE 

SMA 

ECC INC 

PERI 

NCDE 




ITERATICN I 







CONTROL PARAMETERS 
ALPHA 

-56. 3595295 AAO 

BETA 

34.9476030252 


DELTA 

8.09468324C91 

THDOT 

2.407462446428E-02 

TBURN 

2371.64763119 

TASTART 

-61.7465949683 

LAGRANGE MULTIPLIERS 
9A.2ei5A2AC09 

517767.361296 


C. 

182.143038919 

-3404.67365675 

0. 

ORBIT OBTAINED FPCM 
SMA 20455. COC 

PRESENT CONTROLS 
ECC .77662937 

INC 

35.281270 

PER 130.38000 NOD 

303.27000 TAN 

44.620900 

TARGET PARAMETERS 
SMA 

20455. CCCCOOO 

ECC 

.776629372365 


INC 

35.2812702458 

PERI 

130.380000000 

NODE 

303.270CC0CC0 

TAN 

44.6208995261 

DELTA V 1.2172220 







ERRORS IN TARGET VARIABLES 

SMA ECC 

3.2 596290U154E-09 3.9C798 5046681E- 

14 

INC 

-.281270245830 

PERI 

-2.7284041O5319E-12 

NODE 

1.8189894C3546E-12 

TAN 

0. 

CORRECTICNS TO CONTROL VARIABLES 

ALPHA BETA 

4. 53198691813 4. 15337C133L7 


DELTA 

-3.68688135516 

THOOT 

-3.248484222292E-05 

TBURN 

-1.088411549464E-C9 

TASTART 

8.109117220136E-03 

EIGENVALUES OF SECCND PARTIALS OF AUGMENTED 
2.50 7179 829CD+ 13 9.75966656340+06 

FUNCTION 

6.965C759127C+C6 

2.5602547486D+06 

2.49473621C80+C5 

2.23577466720+Cl 


cn 



Oi 


TABLE IV - Continued 


ITERATICN 2 
CCNTPOL PARAMETERS 


ALPHA 

-53.8275426259 

BETA 

39.1C09131564 

DELTA 

5.2C78C186575 

THDOT 

2.404213962206E-02 

TBURN 

2371.64763119 

TASTART 

-61.7384858511 

LAGRANGE MLLTIFLIERS 
93.8639C59660 

558372. 117C21 

24012.208C778 

138.015147139 

-129.568496414 

0. 


ORBIT OBTAINED FFCM PRESENT CONTROLS 

SHA 20A82.25T ECC .77689615 INC 35.035507 PER 130.5497^ NOD 303.16178 TAN A4.570231 


TARGET PARAMETERS 

SMA ECC 

20 A82. 2 51 1873 . 7768961 5C604 

DELTA V 1.211222C 

ERRORS IN TARGET VARIABLES 

SMA ECC 

-27.2571612891 -2. 661782366716E-04 

CORRECTIONS TO CONTROL VARIABLES 

ALPHA BETA 

-.A3587827C093 6 . A7 A 16A7 I 3686E-02 

EIGENVALIES OF SECOND PARTIALS OF AUGMENTED 
2.51513785350+13 9.68C720A6A60+06 


INC PERI 

35.035507C2C7 130. 5 A9736793 


INC PERI 

-3.550702073790E-C2 -.169736792887 


DELTA THDOT 

.362383A971A9 7.086660695794E-05 

FUNCTION 

7.215C75A7430+C6 2.33752503090+06 


NODE TAN 

303.161719A28 44.5702307306 


NODE TAN 

.1C8220572205 0. 


TBURN TASTART 

.6924509229A5 4.973723280717E-03 


3.64442207120+05 2. 22197394 160+Cl 


ITERATION 6 
CONTROL PARAMETERS 

ALPHA BETA DELTA TFDOT TBURN TASTART 

-5A. 2 6156 2A406 39.1678744022 5 . 5696 535 ICC 1 2 . A1 12 3908 17 C3E-02 2372. 370 645A5 -61.7341841275 

LAGRANGE MLLTIPLIERS 

93.8681130452 552C67. 961589 24 3 7 5. 5C9C 6C 9 20.4146420811 -24 .70C 8 1 7726 2 0. 

ORBIT OBTAINED FROM PRESENT CONTROLS 

SMA 20455.00C ECC .77662937 INC 35.COOOOO PER 130.38000 NOD 303.2700C TAN 44.674292 


TARGET PARAMETERS 
SMA 

20455. OCCCOOO 


ECC 

.176629372365 


INC 

35.COOCOOOOOO 


PERI NODE TAN 

130.300000000 303. 2700 COCCO 44.6742915574 


DELTA V 1.2116672 


ERRORS IN TARGET VARIABLES 

SMA ECC INC PERI 

1.979C60471058E-09 1 . 7 76356839400E- 14 9. 09 4947C 1 172 9E- 13 -9. 09494701 7729E-1 3 0. 


NODE 


TAN 

0 . 


CORRECTIONS TO CCNTRCL VARIABLES 

ALPHA BETA 

-4.876S283C2904E-12 -1. 005350177936E -1 I 


DELTA THDOT TBURN TASTART 

2.465918164164E-11 -3 . 87407463049 0E- 15 -4 . 5888 389 2485 8E- 1 1 -1 . 7395239, 23412E-1 1 


EIGENVALUES OF SECOND PARTIALS OF AUGMENTED 
2. 5C 92 86 81 270+ 13 9.6624449378D+06 


FUNCTION 

7.28327668140+C6 


2.34675544790+06 3.28895083710+05 2. 21105692C9D+C1 



TABLE IV - Concluded 


^EHTCN P^PHSC^ ►^TPIX 
1.083E+C7 -3.127E+C6 

1.014E+C7 

1.291E+10 

6.567E+02 -7.C90E+06 

-6.7536+02 

-1. 4376-02 

2.935E+0C 

-2.378E+01 

-3.137F+06 

4.388E+06 

-4.7C9E+C4 

-4.223E+C7 

-7.467E+02 -8.3C1E+04 

1.6366+03 

1.752E-02 

-6.680E+0C 

-6.889E+00 

1.01AE+C7 

-4.7C9E+C4 

1.303E+C7 

1.6616410 

3.448E+C2 -9.257E+C6 

1.9336+02 

-6.C76E-C3 

-5.7CIE-C1 

-4.1676+01 

1.291E+10 

-4.223E+07 

1.661E+10 

2.5C9E+13 

2.416E+C5 -1.194E+10 

1.7316+05 

1. 9276+01 

-1.058E+C3 

-5.2406+04 

6.567E+02 

-7.467E+02 

3. 4486+02 

2.416E+C5 

2.3C5E+C1 2.266E+C3 

-4.7426+01 

-5.C61E-04 

-5.529F-C4 

-1.7316-02 

-7.090E+C6 

-8.3C1E+04 

-9.257E+C6 

-1.1946+10 

2.266E+03 1.469E+07 

1.286E+02 

-1.022E-03 

-4.468E-C1 

-3.115E+01 

-6.753E+C2 

1.636E+03 

1.933E+C2 

1.731E4C5 

-4.142E+C1 1.286E+C2 

0. 

0. 

0. 

0. 

-1.437E-C2 

1. *1528-02 

-6.C76E-C3 

1.927E+C1 

-5.C61E-04 -1.022E-03 

0. 

0. 

0. 

C. 

2.935E+00 

-6.680E+00 

-5.701E-C1 

-1.C58E+03 

-5.529E-C4 -4.468E-01 

0. 

0. 

0. 

0. 

-2.378E+C1 

-6.669E+00 

-4.167E+01 

-5.240E+C4 

-1.731E-02 -3.115E+C1 

0. 

0. 

0. 

C . 

-1.096E+C1 

1.C66E+01 

1.720E+CO 

1.581E+C3 

8.341E-C5 -1.016E+00 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

C. 

c. c . 

0. 

0. 

0. 

c. 

INVERSE 

2.6810-C9 

1.2660-09 

2.9240-09 

1.3040-12 

-2.1650-17 -8.4390-09 

6.6690-05 

-6.1390+00 

-2.293C-01 

2.4390-03 

1.2860-C9 

6.1720-10 

1.4030-09 

6.2580-13 

-1.0390-17 -4.C49D-09 

5.3620-05 

-4.88 7D+0C 

-2.6480-Cl 

3.6000-03 

2. 9240-09 

l.^C 30-09 

3.188C-C9 

1.4220-12 

-2.3610-17 -9.2020-09 

4.4080-04 

-4.1200+01 

1.8590-Cl 

-9.6050-03 

1.3040-12 

6.258C-13 

1.4220-12 

6.3450-16 

-1.0530-20 -4.1050-12 

-3.7970-07 

3.575D-C2 

-6.651C-C6 

-5.2850-06 

-2.1650-17 

-1.C290-17 

-2.3610-17 

-1.0530-20 

1.7490-25 6. 8150-17 

-1.9780-02 

-1.1640+02 

-5. 137C+0C 

-4.3030-03 

-8.4390-09 

-4.C^9C-09 

-9.202C-C9 

-4.1050-12 

6.8150-17 2.6560-08 

-2.8470-06 

8.1710-Cl 

-1.051C-C3 

-1.3020-02 

6.6690-C5 

5.2620-05 

4.4080-04 

-3.7970-C7 

-1.9780-02 -2.6470-06 

-5.7780-01 

5.3610+04 

3.8490+0C 

-7.5890+00 

-6. 1390+C0 

-4. 1670+00 

-4.1200+Cl 

3.5750-C2 

-1.1640+02 8. 1710-01 

5.3610+04 

-5.0610+09 

-3.234D+C5 

7.C310+05 

-2.293D-C1 

-2.6480-01 

1.8590-01 

-6. 6510-06 

-5.1370+00 -1. 0510-03 

3.8490+00 

-3.2340+05 

-8.542C+04 

9.3980+02 

2.4390-C3 

3.6C0D-03 

-9.6050-C3 

-5.2850-06 

-4.3030-03 -1.3020-02 

-7.5890+00 

7.0310+C5 

9.398D+C2 

-1.8410+03 

-1.4210-01 

-6.9890-02 

1.1060-01 

-7.9720-C6 

5.2C60-C3 -1.C650-C2 

-5.3520+00 

5. 0520+05 

-1.105C+C4 

-1.5210+03 

0. 

0. 

0. 

0. 

0. 0. 

0. 

0. 

0. 

C. 


-1.096E+C1 0 

1.066E+01 0 

1.72CE+00 0 

1.581E+03 0 

8,3AlE-05 0 

-1.016E+C0 0 


0 « 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 « 0 


-U421D-01 0 

-6.9890-02 0 

1.1C60-01 0 

-7. 9720-06 0 

5.2C60-03 0 

-1. 0650-02 0 

-5.3520+06 0 

5.C520+05 0 

-1.1050+04 0 

-1.5210+03 0 

-2.058D+C4 0 

0 . 0 


TIME FOP THIS C^SE 29.522 


•3 



tC*' 

00 


TABLE V 

SAMPLE OUTPUT FOR BACKWARD TARGETING MODE 

[Mode = ^ 


SDAT 




NOPT 


i. 1. i, i» 1, 

It 1, 0, 0, 1, 1, 0, 

Kopr 

- 

5f 3, 4, A, 

1. 

AIN 

= 

0,2S^A7E+Clf 0-0, 

0 . 0 , -0.AS4E + 01t 0. 13CCeE+C3t 0.0, 

CONI 

- 

0. 2CA5 5E+C5 , C. 7 1 C6 2S 372 3 65E +CC , 0.35E+02, C-13038E+03, 

0-3C327E+03, 

GS 


-0-55F+C2, 0.4E+02, 

0-56+01, 0.25E-01, C.2AE+0A, 0.^5E+02, 

GL 

= 

0. 1E4C1, O.lE+01 , 

O.IE + Cl, C. lE+01, Q-IE + Cl, C-IE+Oli 

HP 


0.6E+CCf 0-6E+0C, 

0.6E + 00, C.6E-03, O.lE+01 , 0.6E + 00t 

VI 

= 

0-15E+02, 0.15E+02 

, 0.15E+02, 0.15E-02, C.5E+02, 0.15E+02 

NSTEPS 

= 

8, 


PASS 


0.32C6S5E4CA, 


DMASS 

= 

-0. An2E+CC, 


THR 

= 

0.132A5E+C1, 


PERJD 

= 

It 


REFJO 

= 

0.0, 


PER 


It 


SLAT 

=: 

It 


SLON 


It 


ERR 


O.U-07, 


MAXIT 


50, 


UMARS 

= 

O.^28284E+05t 


MODE 


3t 


BOUND 


0.0, 0.18E403, 



$ENO 



TABLE V - Continued 


VIKING TARGETING 

♦♦♦♦♦BACKWARD TARGETING ( MOD E=3 ) ♦♦♦♦♦ 

INITIAL CCMC 

SMA 20455.00C ECC .7T662937 INC 35.CCCCOO 

THE CCNTFCL VARIABLES ARE ALPHA BETA DELTA 

THE TARGET VARIABLES ARE VINE OECSV RTASV 


ITERATICN 1 
CONTROL PARAMETERS 

ALPHA BETA DELTA 

-55,OCOOQCCOOO 40 . CCCOCOOOOO 5 .CCCCCCOCOCC 

LAGRANGE MLLTIPLIERS 

l.COCOOOCCOOO 0. 0. 

ORBIT OBTAINED FRCM PRESENT CCNTBOLS 
SMA -4854.6847 ECC 1.9907407 INC 35.441364 

TARGET PARAMETERS 

VINE B^R B+T 

2.97C198C7011 -48 C4. OC 1 83 1 78 6837.73257053 

DELTA V 1.2355143 

ERRORS IN TARGET VARIABLES 

VINE B^R B^T 

-2. 5498070 1C624E-02 0. 0. 

CORRECTICNS TQ CONTROL VARIABLES 

ALPHA BETA DELTA 

.462710CCZ179 .457260657010 -2 . C 33737 7455 0 

FINAL MASS 2061.67C00C 

ETGENVALLES OF SECOND PARTIALS OF AUGMENTED FUNCTION 

8.33241465 41D+01 2. 40 839 7 34 5 7D<- 00 1 . 9 99 9 94 1 81 6D + 00 


ANALYSIS PROGRAM (VITAP) 

PER 130.38000 NOD 303.27000 TAN 45.0CQ000 

THDCT TBLRN TASTART 


THDOT TBURN TASTART 

2.5COCOOOOOOOOE-02 2400. OOCCOOOO 45.0000000000 


l.COOCOOOOOOO UOOOOCOOOCCO 0. 


PER 129.36505 NOD 302.82894 TAN 297.86829 


DECSV RTASV TAN 

354.674089349 130.355024490 297.868291354 


OECSV RTASV TAN 

.385910651461 -.275024490129 0. 


THDOT TBURN TASTART 

L.500000000000E-03 -4.40296290313 -.155157459117 


-3. 72 36 3 23400D4-00 -6 . 13476980 160+00 -6 . 279 17766 2 5D+06 


CD 



CJl 

o 


TABLE V ~ Continued 


ITERATICN 2 


CCNTRCL PARAMETERS 
ALPHA 

-54.5172899970 

BETA 

4C.45726C6570 

DELTA 

6.66973 108272 

LAGRANGE MULTIPLIERS 
2.305861^7833 

0. 

C. 

ORBIT OBTAINED FPCM 
SMA -4861.5682 

PRESENT CCNTPCLS 
ECC 1.9850033 

INC 35.456375 


THDOT 

TBURN 


TASTART 

2.3500CCCOOOOOB-02 

2395.597C3710 

45. 

1551574591 

.979964684660 

.990405855121 

0. 


129.41727 NOD 

302.83375 

TAN 297 

.86383 


TARGET PARAMETERS 
VINF 

2.96809456648 

B*R 

-4796.90561222 

04 T 

6823.25774489 

DELTA V 1.23Z6658 



ERRORS IN TARGET VARIABLES 

VINF B*R 

-2.339456647724E-C2 0. 

B + T 

0. 

CORRECTICNS TO CCNTRCL 
ALPHA 

.572259724437 

VARIABLES 

BETA 

.545786574855 

DELTA 

-2.05866091386 


FINAL MASS 2C63.771C9A 


DECSV RTASV TAN 

354.6 897095C2 130.333321173 297.863830<*83 


DECSV RTASV TAN 

.370290498136 - . 25333 1 1 72693 0. 


THDOT TBUR^ TASTART 

1.500CC0OOO000E-03 -3.16638921221 -.116688788908 


EIGENVALLES OF SECOND PAPTIALS OF AUGMENTED FUNCTION 

5.C960371E67D+01 2. 000C053276D+OO 1 . 5 2 33947 392C ♦ CO -4.79416468630+00 -7. 79 154682C80+00 - 1 . C 1 2677 57 C6D+C7 


ITERATICN 46 


CONTRCL PARAMETERS 

ALPHA 

BETA 

DELTA 

-58.97538(4989 

34.5884174529 

16.291160' 

LAGRANGE MULTIPLIERS 

-7460119.32293 

0. 

C . 

ORBIT OBTAINED FPCM 
SMA -4939.1218 

PRESENT CCNTPOLS 
ECC 1.9596038 

INC 35.C5C806 


THDOT TBURN TASTART 

1.9C1992734839E-02 2367.96688861 44.5623489078 

9783.80393437 8895.53598726 0. 

PER 129.3C790 NOO 303.00276 TAN 297.70925 


TARGET PARAMETERS 
VENF 

2. 94469999956 


B^R 

-4743.85826518 


B*T 

6839.48432114 


DECSV 

355.060000003 


RTASV 

130.079999996 


TAN 

297.709246057 


DELTA V 


1.2HE561 


ERRORS IN TARGET VARIABLES 
VINF 

4.429C97C89C87E-1C 0. 


B*R 


B4T 


0 . 


DECSV 

-3.044988261536E-09 


RTASV 

4.349203663678E-09 


TAN 


0 . 


CORRECTICNS TO CCNTRCL VARIABLES 

ALPHA BETA 

1.4 7 73480«e775E-C9 1 . 1 * 309 70 04863E-C9 


DELTA 

3.4CE362694156E-C9 


THDOT 

1.168156685974E-12 


TBURN TASTART 

6.957999209935E-07 -9.931706836794E-09 


FINAL MASS 


2C76. 956201 


EIGENVALUES OF SECOND PARTIALS OF AUGMENTED FUNCTION 

1.99E43231‘2D + 13 1. 0409CJ 748 80 + 07 5. 8 19 7 7255 02D + C6 1.92 532 6 58810+06 2 . 364 30 1 1 5 20D + O 5 3.788190 229 20+00 



NASA-Langley, 1971 30 L-7440 


TABLE V - Concluded 


KEWTCN R^PHSON P^TRIX 


9.6A2E+C6 

-3.036E+06 

9.00AE+C6 

9.753E+C9 

-5.326E+01 

5.462F+06 

1.340E-02 

0. 

0. 

-5.643E+00 

1.7456+01 

0. 

-3.036E+C6 

5.264E+06 

-3.51CE+04 

-3.490E+C7 

-U239E+01 

1.989E+C4 

-l.lllE-02 

0. 

0. 

-5.290E+00 

-3.496E + C0 

0. 

9.004E4C6 

-3.510E+0^ 

1.C91E+C7 

l.ie3E + lC 

-7.768E+C1 

6.648E+C6 

5.376E-03 

0. 

0. 

-1.201E+01 

1.7726+01 

0. 

9.753E4C9 

-3.^90E+07 

1.183E+10 

1.998E413 

-6.205E+C5 

6.970E+09 

5.571E+00 

0. 

0. 

-1.218E+04 

1.8C6E+04 

0. 

-5.326E+CI 

-1. 239E+C1 

-7.760E+C1 

-6.205E+C5 

3.918E+00 

4.6C8E+02 

6.352E-04 

0. 

0. 

-1.4806-03 

1.961E-03 

0. 

5.^62E+06 

1.989E+0A 

6.6A8E+06 

6.970E+C9 

4.6C8E4C2 

6.7606+06 

-2.079E-03 

0. 

0. 

4.6176+00 

-6.822E+00 

0. 

1.340E-C2 

-1. 111E-C2 

5.376E-C3 

5.571E4C0 

6. 352E-C4 

-2.C79E-03 

0. 

0. 

0. 

C. 

0. 

0. 

a. 

3. 

0. 

C. 

C. 

C. 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

0. 

0. 

0. 

C. 

C. 

0. 

0. 

0. 

c. 

0. 

0. 

-5.6A3E+00 

-5.2SCE+C0 

-1.201E+C1 

~l.218E+C^ 

-1.48CE-03 

4.617E+00 

0. 

0. 

0. 

c. 

0. 

0. 

1.745E+01 

-3. AS6E+00 

1.772E+01 

1.8C6E4C^ 

1.96 1E-C3 

-6.622E+C0 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

0. 

0. 

0. 

C. 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

INVERSE 

1. 268D-C6 

1.C23D-C6 

♦l.OAAD-06 

-1.521D-12 

-4.9210-15 

2.C70C-09 

4.1340-02 

0. 

c. 

2.C52D-CI 

1.400D-01 

0. 

1.023D-C6 

8.262D-C7 

-8.^A1C-C7 

3.685D-IA 

-3.9260-15 

1. 5550-09 

3.5460-03 

0. 

0. 

3.5150-02 

2.4600-02 

0. 

-1.04^0-06 

-8.A^1D-C7 

i.oieo -06 

-I.A20D-10 

-2.951C-16 

3.1930-00 

-3.1250-01 

c. 

c. 

-1.8850-01 

-8.6810-02 

0. 

-1.521C-12 

3. tE5C-lA 

-1.A20C-10 

1.389D-13 

5.4170-18 

-4.9580-12 

1.3030-04 

0. 

0. 

1.6530-C6 

-1.529C-C6 

0. 

-^.9210-15 

“3.926C-15 

-2.951C-16 

5.^170-18 

3.6870-22 

3.C060-15 

1. 5750+03 

0. 

0. 

-2.066D+00 

-1.8780+00 

0. 

2.07C0-C9 

1.555D-09 

3.193C-C8 

-4.9580-12 

3.0C60-15 

7.4300-08 

9.0200-02 

0. 

0. 

2.0860-02 

-3.116D-02 

0. 

4.134D-C2 

3.5A60-03 

-3.1250-Cl 

1.3030-C4 

1.57504C3 

9.0200-02 

-9.0820+06 

c. 

0. 

-2.8720+04 

4.1070+C4 

0. 

0. 

0. 

0. 

0. 

C. 

C. 

0. 

c. 

c . 

C « 

0. 

0. 

0. 

0. 

C. 

0. 

c. 

0. 

0. 

0. 

0. 

c. 

0 . 

0. 

2.052D-01 

3. 5 15C-02 

-1.885C-01 

1.6530-06 

-2.0660+00 

2.C86D-02 

-2.8720+04 

0. 

0. 

-5.7570+04 

-3. 6060+04 

0. 

1.4C0D-C1 

2.^600-02 

-8.681C-C2 

-1.529D-C6 

-1.8760+00 

-3.116C-C2 

4.1070+04 

0. 

0. 

-3.6C60+04 

-2.9400+C4 

0. 

0. 

0. 

0. 

C. 

0. 

C. 

0. 

0. 

0. 

C . 

0. 

0. 
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